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Neutron stars are composed of the densest form of matter known 
to exist in our universe, and thus provide a unique laboratory for 
exploring the properties of cold matter at supranuclear density. 
Measurements of the masses or radii of these objects can strongly 
constrain the neutron-star matter equation of state, and conse- 
quently the interior composition of neutron stars [TJ, [2]. Neutron 
stars that are visible as millisecond radio pulsars are especially 
useful in this respect, as timing observations of the radio pulses 
provide an extremely precise probe of both the pulsar's motion and 
the surrounding space-time metric. In particular, for a pulsar in a 
binary system, detection of the general relativistic Shapiro delay 
allows us to infer the masses of both the neutron star and its binary 
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companion to high precision [3, 4|. Here we present radio timing ob- 
servations of the binary millisecond pulsar PSR J1614— 2230, which 
show a strong Shapiro delay signature. The implied pulsar mass 
of 1.97±0.04 M Q is by far the highest yet measured with such cer- 
tainty, and effectively rules out the presence of hyperons, bosons, 
or free quarks at densities comparable to the nuclear saturation 
density. 

In the accepted "lighthouse model" description of radio pulsars, a rapidly 
spinning neutron star (NS) with a strong magnetic field (10 8 -10 15 G) emits a 
beam of radiation that is typically misaligned with the spin axis. A broad- 
band, polarized pulse of radio emission is observed once per rotation if this 
beam crosses the Earth-pulsar line of sight. The extraordinary rotational sta- 
bility of pulsars permits the precise measurement of a number of systematic 
effects that alter the arrival times of the radio pulses at Earth, a procedure 
referred to as pulsar timing. In the case of binary millisecond pulsars, which 
are the most stable pulsars with orbital companions, even typically subtle 
effects such as the general relativistic Shapiro delay can be revealed by tim- 
ing. The Shapiro delay is an increase in light travel time through the curved 
space-time near a massive body. In binary pulsar systems that have highly 
inclined (nearly edge-on) orbits, excess delay in the pulse arrival times can 
be observed when the pulsar is situated nearly behind the companion dur- 
ing orbital conjunction. As described by general relativity, the two physical 
parameters that characterize the Shapiro delay are the companion mass and 
inclination angle. In combination with the observed Keplerian mass function, 
the Shapiro delay offers one of the most precise methods to directly infer the 
mass of the NS. In turn, any precise NS mass measurement limits the equa- 
tions of state (EOS) available to describe matter at supranuclear densities. 
The discovery of a NS with mass significantly higher than the typical value 
of ~1.4M would have a major impact on the allowed NS EOS as well as 
additional implications for a wide range of astrophysical phenomena [5]. 

PSR J1614— 2230 was originally discovered in a radio survey of uniden- 
tified EGRET gamma-ray sources using the Parkes radio telescope [6J. The 
spin period P is 3.15 ms, and initial timing with Parkes showed the pulsar to 
be in a binary system with an 8.7-day orbital period and a companion of mass 
M2 > 0.4 M . The system was noted as having a higher companion mass 
than is typical for fully-recycled (P < 10 ms) pulsars, which predominantly 
have helium white dwarf (WD) companions with masses of ~0.1— 0.2 M . 
Furthermore, the orbital period is shorter than expected given the massive 
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companion [7J . These facts hinted at a possible non-standard evolutionary 
history for the binary system (see below). 

Following the inital discovery observations, J1614— 2230 was observed 
regularly as a test source for several observing projects at the NRAO Green 
Bank Telescope (GBT)[8j3. These data provide a continuous long-term tim- 
ing record from mid-2002 to the present, and they show a marginally signifi- 
cant Shapiro delay signal. However, the data quality and arbitrary scheduling 
with respect to the binary's orbital phase made it impossible to obtain mean- 
ingful mass and inclination measurements from those data alone. In March 
2010, we performed a dense set of observations of J1614— 2230, timed to fol- 
low the system through one complete orbit with special attention paid to the 
orbital conjunction, where the Shapiro delay signal is strongest. These data 
were taken with the newly built Green Bank Ultimate Pulsar Processing 
Instrument [9] . GUPPI coherently removes interstellar dispersive smearing 
from the pulsar signal, and integrates the data modulo the current appar- 
ent pulse period, producing a set of average pulse profiles, or flux-versus- 
rotational-phase lightcurves. We observed an 800 MHz wide band centered 
at a radio frequency of 1.5 GHz. The raw profiles were then polarization- 
and flux-calibrated and averaged into 100-MHz, 7.5-minute intervals using 
the PSRCHIVE software package [TO]. From these, pulse times of arrival were 
determined using standard procedures, with a typical uncertainty of ~1 fis. 

The measured arrival times are used to determine key physical parameters 
about the neutron star and its binary system by fitting them to a compre- 
hensive timing model which accounts for every rotation of the neutron star 
over the time spanned by the fit. The model predicts at what times pulses 
should arrive at Earth, taking into account pulsar rotation and spin-down, 
astrometric terms (sky position, proper motion), binary orbital parameters, 
interstellar dispersion, and general relativistic effects such as the Shapiro de- 
lay (see Table [1]). The observed arrival times are compared with the model 
prediction, and best-fit parameters are obtained via x 2 minimization, using 
the TEMP02 software package [11 J*[ The post-fit residuals, i.e. the differ- 
ence between the observed and model-predicted pulse arrival times, effec- 
tively measure how well the timing model describes the data, and are shown 
in Figure [H We included both the older long-term datasets and our new 

'''The National Radio Astronomy Observatory is a facility of the U.S. National Science 
Foundation (NSF), operated under cooperative agreement by Associated Universities, Inc. 
*We also obtained consistent results using the original TEMPO package. 
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GUPPI data in a single fit. The long-term data determine model parameters 
with characteristic timescales longer than a few weeks (e.g. spin-down rate 
and astrometry), while the new data best constrain parameters on timescales 
of the orbital period or less. Both are necessary to reduce covariance between 
particular model parameters. 

In addition to the physical timing model parameters listed in Table [TJ 
the fit included arbitrary time offsets ("jumps") between different observing 
systems, and allowance for a time-variable interstellar dispersion measure 
(DM), a quantity proportional to the total free electron column density along 
the line of sight. A single DM value was fit for the long-term data set. The 
newer, more precise GUPPI data required a separate DM be fit for each day. 
DM is expected to vary on 1-day and longer timescales due to the relative 
motions of the pulsar and Earth through the interstellar medium |12j. We 
also investigated the possibility of DM variation on minute to hour timescales 
which could indicate the presence of excess gas in the binary system itself. If 
present, especially near orbital conjunction, this effect could bias the Shapiro 
delay measurement. However, we found no evidence of this in our data - 
we limit any such orbital DM variation to a level of < 2 x 10~ 4 pccm -3 , 
corresponding to delays of <0.4/xs at a frequency of 1.5 GHz. Furthermore, 
we observed orbital conjunctions at three different epochs over a span of 2 
months, and all show consistent timing results without needing to model fast 
DM variation. 

As shown in Figure [TJ the Shapiro delay was detected in our data with 
extremely high significance, and must be included to correctly model the 
arrival times of the radio pulses. However, estimating parameter values and 
uncertainties can be tricky due to the high covariance between many orbital 
timing model terms, especially the Shapiro-derived companion mass (M 2 ) 
and inclination angle (i) [T5] . In order to obtain robust error estimates, we 
used a Markov chain Monte Carlo (MCMC) approach to explore the post-fit 
X 2 space and derive posterior probability distributions for the model param- 
eters (see Figure [2]). In contrast with the common method of mapping x 2 
projections in a reduced-dimensional parameter space[HJ G2, MCMC is able 
to efficiently explore all fit dimensions simultaneously. Our final results for 
all timing model parameters, with MCMC error estimates, are given in Ta- 
ble [TJ Due to the high significance of the detection, our MCMC procedure 
and the standard x 2 fit produce similar uncertainties. 

From the detected Shapiro delay we measure a companion mass of 0.500 ± 
0.006 M , which implies that the companion is a helium-carbon-oxygen white 
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dwarf, as helium WD stars have masses of < O.45M [T2]. Considering the 
8.7-day orbital period, a likely formation mechanism for systems similar to 
J1614— 2230 has been proposed [T6l IT7] where the original companion star is 
~3M . J1614— 2230 is peculiar, however, in being fully recycled unlike the 
mildly-recycled pulsars (P ~ 10 — 30 ms) found in similar binary systems. 

The Shapiro delay also shows the binary system to be amazingly edge-on, 
with an inclination of 89.17° ± 0.02°. This is the most inclined pulsar binary 
system currently known. The amplitude and sharpness of the Shapiro delay 
increase rapidly with increasing binary inclination and the overall scaling of 
the signal is linearly proportional to the mass of the companion star. Thus 
the unique combination of the high orbital inclination and massive WD com- 
panion in J1614— 2230 cause a Shapiro delay amplitude orders of magnitude 
larger than for most other millisecond pulsars with 0.1— 0.2 M companions 
and inclinations <75°. In addition, the excellent timing precision achievable 
from the pulsar with the GBT and GUPPI provide a very high signal-to-noise 
ratio measurement of both Shapiro delay parameters within a single orbit. 

The standard Keplerian orbital parameters, combined with known com- 
panion mass and orbital inclination, fully describe the dynamics of a "clean" 
binary systenj§ under general relativity and therefore also determine the pul- 
sar's mass. We measure a NS mass of 1.97 ± 0.04 M , by far the highest pre- 
cisely measured NS mass determined to date. In constrast with X-ray based 
NS mass/radius measurements [T8~| . the Shapiro delay provides no informa- 
tion about the NS radius. However, unlike the X-ray methods, our result is 
nearly model- independent, as it depends only on general relativity being an 
adequate description of gravity. In addition, unlike statistical pulsar mass de- 
terminations based on measurement of the advance of periastron[T9l 1201 [21] . 
pure Shapiro delay mass measurements involve no assumptions about classi- 
cal contributions to periastron advance or the distribution of orbital inclina- 
tions. 

The mass measurement alone of a 1.97 M NS significantly constrains the 
NS EOS, as shown in Figure [3j Any proposed EOS whose mass-radius track 
does not intersect the J1614— 2230 mass line is ruled out by this measurement. 
The EOSs which produce the lowest maximum masses tend to be those which 
predict significant softening past a certain central density. This is a common 
feature of models that include the appearance of "exotic" hadronic matter 
such as hyperons or kaon condensates at densities of a few times the nuclear 

§A system comprised of two stable compact objects. 
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saturation density (n s ), for example GS1 and GM3 in FigureEJ All such EOSs 
are ruled out by our results. Our mass measurement does not completely rule 
out condensed quark matter as a component of the NS interior [22J, but it does 
constrain quark matter model parameters, and excludes some strange quark 
matter models[5]. For the range of allowed EOS lines presented in Figure [31 
typical values for the physical parameters of J1614— 2230 are a central baryon 
density between 2— 5n s and a radius between 11 — 15km, only 2—3 times 
the Schwarzschild radius for a 1.97 M Q star. Based on EOS-independent 
analytic solutions of Einstein's equations, our mass measurement also sets 
an upper limit on the maximum possible mass density of cold matter [23]. of 
< 4 x 10 15 g cm" 3 ~10 n s . 

The proposed formation mechanism for J1614— 2230-like systems [TBI H77] 
predicts that the neutron star would only accrete a few hundredths of a solar 
mass of material, much less than the ~O.6M needed to bring a neutron 
star born at 1.4 M up to the observed pulsar mass. This implies either that 
the transfer of large amounts of mass onto the neutron star from the sec- 
ondary is possible or alternatively that some neutron stars are born massive 
(~1.9M ). As the Shapiro delay has not been detected in most millisecond 
pulsar systems, our mass measurement for J1614— 2230 suggests that many 
of these other systems may also harbor NSs with masses well above 1.4 M . 
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Figure 1: Shapiro delay measurement for PSR J1614— 2230. Each panel 
shows timing residual - excess delay not accounted for by the timing model 
- as a function of the pulsar's orbital phase. The top panel shows the full 
magnitude of the Shapiro delay if it is not included in the timing model 
but with all other model parameters fixed at their best-fit values. The solid 
line shows the Shapiro delay functional form, and the red points are timing 
measurements from our GBT/GUPPI data set. The diagrams inset in this 
panel show a top-down schematic view of the binary system at orbital phases 
0.25, 0.5, and 0.75. The neutron star is shown in red, the white dwarf 
companion in blue, and the emitted radio beam (yellow) points towards the 
Earth. At orbital phase 0.25, the Earth-pulsar line of sight passes nearest 
to the companion (~240,000 km), producing the sharp peak in pulse delay. 
We found no evidence for any kind of pulse intensity variations, as from 
an eclipse, near conjunction. The middle panel shows the best-fit residuals 
obtained using an orbital model that does not account for general-relativistic 
effects. That the residuals deviate significantly from a random, Gaussian 
distribution of zero mean shows that the Shapiro delay must be included to 
properly model the pulse arrival times, especially at conjunction. In addition 
to the red GBT/GUPPI points, the gray points show the previous "long- 
term" data set described in the text. The dramatic improvement in data 
quality is apparent. Finally, the bottom panel shows the post-fit residuals 
for the fully relativistic timing model (including Shapiro delay), which have 
a root-mean-square residual of 1.1 /is and reduced x 2 value of 1.4 with 2165 
degrees of freedom. 
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Figure 2: Results of the Markov chain Monte Carlo (MCMC) error analysis. 
The left panel grayscale shows the 2-D posterior probability density func- 
tion (PDF) in the M2-1 (companion mass versus orbital inclination) plane 
computed from a histogram of MCMC trial values. The ellipses show 1- and 
3-cr contours based on a Gaussian approximation to the MCMC results. The 
right panel shows the PDF for pulsar mass derived from the MCMC trials. 
In both cases the results are very well described by normal distributions due 
to the extremely high signal-to-noise ratio of our Shapiro delay detection. 
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Figure 3: Neutron star (NS) mass-radius diagram. The plot shows non- 
rotating mass versus physical radius for several typical NS equations of state 
(EOS) [25] . The horizontal bands show the observational constraint from our 
J1614— 2230 mass measurement of 1.97±0.04 M Q , similar measurements for 
two other millsecond pulsars [21 EE], and the range of observed masses for 
double NS binaries[2]. Any EOS line that does not intersect the J1614— 2230 
band is ruled out by this measurement. In particular, most EOS curves in- 
volving exotic matter, such as kaon condensates or hyperons, tend to predict 
maximum NS masses well below 2.0 M Q , and are therefore ruled out. 
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Ecliptic Longitude, A (deg) 245.78827556(5) 

Ecliptic Latitude, (3 (deg) -1.256744(2) 

Proper Motion in A (mas yr" 1 ) 9.79(7) 

Proper Motion in f3 (mas yr -1 ) -30(3) 

Parallax (mas) 0.5(6) 

Pulsar Spin Period (ms) 3.1508076534271(6) 

Period Derivative (s/s) 9.6216(9) x 10" 21 

Reference Epoch (MJD) 53600 

Dispersion Measure (pc cm -3 ) 34.4865* 

Orbital Period (days) 8.6866194196(2) 

Projected Semi-Major Axis (lt-s) 11.2911975(2) 

1 st Laplace Parameter, esinu; 1.1(3) xl0~ 7 

2 nd Laplace Parameter, e cos u -1.29(3) xl0~ 6 

Companion Mass (M ) 0.500(6) 

Sine of Inclination Angle 0.999894(5) 

Epoch of Ascending Node (MJD) 52331.1701098(3) 

Span of Timing Data (MJD) 52469-55330 

Number of TOAs 2,206 (454 / 1,752)* 

RMS TOA Residual (jjs) 1.1 

Right Ascension (J2000) 16 h 14 m 36*5051(5) 

Declination (J2000) -22° 30' 31'.'081(7) 

Orbital Eccentricity, e 1.30(4) xl0~ 6 

Inclination Angle (deg) 89.17(2) 

Pulsar Mass (M ) 1.97(4) 

DM-derived Distance (kpc) 1.2 § 

Parallax Distance (kpc) >0.9 

Surface Magnetic Field (10 8 G) 1.8 

Characteristic Age (Gyr) 5.2 

Spin-down Luminosity (10 34 ergs _1 ) .... 1.2 

Average Flux Density at 1.4 GHz (mJy) 1.2* 

Spectral Index, 1.1-1.9 GHz -1.9(1) 

Rotation Measure (rad m~ 2 ) -28.0(3) 



Table 1: Physical parameters for PSR J1614— 2230: Timing model parame- 
ters (top); quantities derived from timing model parameter values (middle); 
radio spectral and interstellar medium properties (bottom). Timing values 
in parentheses represent l-o~ uncertainty in the final digit, as determined by 
Markov chain Monte Carlo error analysis. t : These quantities vary stochas- 
tically on >l-day timescales. Values presented here are the average for our 
GUPPI data set. p. Shown in parentheses are separate values for the long- 
term (first) and new (second) data sets described in the text. §: Calculated 
using the NE2001 pulsar distance mo del [24"]. 



